Variation of Periods Modulo p in Arithmetic Dynamics 
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Abstract. Let ip : V ^ V he d. self-morphism of a quasiprojec- 
tive variety defined over a number field K and let P e V{K) be 
a point with infinite orbit under iteration of For each prime p 
of good reduction, let mp{(p,P) be the size of the (/?-orbit of the 
reduction of P modulo p. Fix any e > 0. We show that for almost 
all primes p in the sense of analytic density, the orbit size mp [ip, P) 
is larger than (log N i 



Introduction 



Let 



be a morphism of degree d defined over Q and let P G P^(Q) be a 
point with infinite forward orbit 

0^(P) = {P,^(P),Vp2(P),...}. 

For all but finitely many primes p, we can reduce ip to obtain a mor- 
phism 

: P^^ — . P^^ 

whose degree is still d. We write mp{ip,P) for the size of the orbit of 
the reduced point P = P mod p, 

m,(^,P) = #0^^(P). 

(For the remaining primes we define mp{(p, P) to be oo.) 

Using an elementary height argument (see Corollary 10), one can 
show that 

mp(y9, P) > (i log logp + 0(1) for all p, 

but this is a very weak lower bound for the size of the mod p orbits. 
Our principal results say that for most primes p, we can do (almost) 
exponentially better. In the following result, we write S{V) for the 
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logarithmic analytic density of a set of primes V. (See Section 1 for 
the precise definition of 5 and the associated lower density 5_.) 

Theorem 1. With notation as above, we have the following: 

(a) For all^ <1, 

6{p : mp{ip,P) > (logp)^} = 1. 

(b) There is a constant C — C{N, ip, P) so that for all e > 0, 

S\^p : mp{<^, P) > e logp} > 1 — Ce. 

More generally, we prove analogous results for any self-morphism 
(f : V ^ V of a quasiprojective variety V defined over a number 
field K. See Section 1 for the basic setup and Theorem 3 for the 
precise statement of our main result. 

The proof of Theorem 1, and its generalization Theorem 3, proceeds 
in two steps. In the first step we prove that there is an integer D{m) 
satisfying loglogD(m) <S m with the property that 

mp{(f, P) <m if and only if p\D{m). 

This is done using a height estimate for rational maps (Proposition 4) 
and a height estimate for arithmetic distances (Proposition 6). The 
second part of the proof uses the first part to prove an analytic es- 
timate (Theorem 11) of the following form: for all A > 1 there is a 
constant C — C{(p, A) so that 

for all. >0. (1) 



pQsmp{ip,P)^ — gl/\ 

p prime 

The inequality (1) is a dynamical analogue of the results in [9], which 
treated the case of periods modulo p of points in algebraic groups. 

Theorem 1 says that for most p, the mod p orbit of P has size 
(almost) as large as log p. If were a random map, we would expect 
most orbits to have size on the order of ^#P^(Fp) ^ p^/^. In Section 6 
we do some experiments using quadratic polynomials (pc{z) = + c 
which seem to suggest that mp{(pc, «) grows slightly slower than y/p. 

1. Notation and Statement of Main Result 

In this section we set notation for our basic objects of study, give 
some basic definitions, and state our main result. We start with the 
dynamical setup. 

Definition. Let K/Q be a number field, let V C be a quasiprojec- 
tive variety defined over K, and let 

ip:V — >V 
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be a morphism defined over K. Fix a point P G V{K) whose forward (p- 
orbit 

is infinite. For each prime ideal p of X at which V and (p have good 
reduction, let 

rrip — mp{(f, P) — size of the (.p-orbit of P in V^(Fp). 

If V has bad reduction at p, we set rup — oo. 

Remark 2. The question of good versus bad reduction requires choosing 
models for V and ip over the ring of integers of K. However, different 
choices affect only finitely many of the 'mp{(p,,P) values, and thus 
have no effect on our density results. We thus assume throughout that 
a particular model has been fixed. 

We next define the analytic density that will be used in the statement 
of our main result. 

Definition. Let K/Q be a number field with ring of integers Rk- For 
any set of primes V C Spec(i?i<:), define the partial ^-function for V by 



p<=V 

This Euler product defines an analytic function on Re(s) > 1. As usual, 
we write Cxis) for the ^-function of the field K. Then the {logarithmic 
analytic) density of V is given by the following limit, assuming that 
the limit exists: 

'"•'^> .-1+ <iiogCA-(.s-) .-1+ c;f(s)/c*-(s) 

Expanding the logarithm before differentiating and using the fact that 
Ck{s) has a simple pole at s = 1, it is easy to check that the density is 
also given by the formula 

We similarly define upper and lower densities S{V) and S{V) by re- 
placing the limit with the limsup or the liminf, respectively. 

With this notation, we can now state our main result. 

Theorem 3. Let K/Q be a number field and let (p : V/K ^ V/K, 
and P e V{K) be as described in this section. Further, let mp{p, P) 
denote the size of the (p-orbit of P in V(¥p). 



4 JOSEPH H. SILVERMAN 

(a) For all '-f < 1 we have 

S{p G Spec Rk : mp{ip,P) > (logNp)^} = 1. 

(b) There is a constant C = C{K, V, ip, P) so that for all e > 0, 

^{p e Spec Rk '■ rrip {^,P) > elogNp} > 1-Ce. 

2. Height and norm estimates 

In this section we prove various estimates for heights and norms that 
will be needed for the proof of our main result. To ease notation, for 
the remainder of this paper wc fix the number field K/Q and write No 
for the K/Q norm of a fractional ideal a oi K. 

Proposition 4. With notation as in Section 1 and Theorem 3, there 
are constants 

d = d{V, ip)>2 and C = C(V, ip)>0 

so that 

h{^p''{Q)) < d!^{h{Q) + C) for alln>0 and all QeV(K). 

Proof. We are given that </? is a morphism on V, but note that V is 
only quasiprojective, i.e., y is a Zariski open subset of a Zariski closed 
subset of P^. We write V as a union of open subsets Vi, . . . ,Vt such 
that on each Vi we can write 

— v\vi — [Fio, Fil, ■ ■ ■ , FiN], 

where the F^j are homogeneous polynomials and such that 

Fio, Fil, ... , FiN do not simultaneously vanish on Vi. 

We may view </?i as a rational map (pi : ¥^ — -»• of degree di — 
deg Fij. Letting Zi C P^ be the locus of indeterminacy for the rational 
map (fii, we have the elementary height estimate 

h{ipi{Q)) < dih{Q) + C((/?i), valid for all Q e P^(i?) \ Zi. (2) 

(See [6, Theorem B.2.5(a)].) By construction, 

y n Zi n ^2 n • • • n = 0, 

so for all Q e V{K) we obtain the inequality 

h{^{Q)) = h(^(pi{Q)) for any i with Q ^ Zi, 

< max dih{Q) + C {(fi) from (2), 

i with Q ^ Zi 

< max dih{Q) + max C{ipi). 

\<i<n l<i<n 
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Setting 

d = max{2,di, . . . ,dt} and C = max{C((^i), . . . , C((^t)}, 
we have 

h{^{Q)) < dh{Q) + C for all Q e V{K). 

Applying this iteratively yields 

h{<p''{Q)) < (ThiQ) + (1 + d + • • • + d"-^)C < d^{h{Q) + C), (3) 

(note that d > 2 by assumption) which is the desired result. □ 

Remark 5. If <^ : ^ is a finite morphism, then in the statement 
of Proposition 4 we can take d — max{deg </?, 2}. More precisely, in 
this situation a standard property of height functions [6, B. 2. 5(b)] gives 
upper and lower bounds, 

/i((/p"(g)) = d"(MQ) + o(i)). 

For maps of degree 1, the middle inequality in (3) yields the stronger 
estimate 

hi^^'iQ)) <h{Q)+Cn. 

Tracing through the proofs in this paper, this would give an exponential 
improvement in our results. For example, consider a linear map Lp{z) = 
az with a e Q*, so mp{ip, 1) is the order of a in the multiplicative 
group F*. Then in place of (1) we would obtain 

y i"gp^ <g + o(i), (4) 

p prime r \ i ' / 

and this would allow us to replace Theorem 1(b) with 

6{p : mp{ip, 1) > > 1 - 2e. (5) 

However, we will not pursue the degree 1 case, because both (4) and (5) 
are special cases of the general results on algebraic groups proven in [9], 
see in particular [9, equation (3)] and the remark following [9, Theo- 
rem 4.2]. 

Proposition 6. Let K/Q be a number field and let 

ao, . . .,aN,Po, ...,Pn & K 

be elements of K with at least one ai and at least one (5i nonzero. 
Define fractional ideals 

21 = (Q;o,...,Q;iv), ^ = (/3o, ■ ■ ■ ,/3iv), £> = (oii/Jj - Q:jA)o<j<j<Ar. 
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Also let A = [ao, . . . , ajv] e F^{K) and B = [(3o, . . . , (3n] G F^{K), 



where the 0{1) depends only on N. {Here h is the absolute logarithm 
height on P^(Q).) 

Proof. We give a proof using the machinery of heights relative to closed 
subschemes developed in [11], although for the specific result that we 
need, some readers may prefer to write down a direct proof. We assume 
that the absolute values on K have been normalized so as to obtain 
the absolute height, i.e., the height of a point P — [xq, . . . , xn] is given 
by h{P) = E^GMx ^^{-v{xi)}. 

Since A ^ B, there are indices i and j such that aif3j ^ o.j(3i. 
Relabeling the coordinates, we may assume without loss of generality 
that 7^ cn/3o- We define subvarieties D and A of x by 

D = {xqUi = XiUq} and A = {xii/j = Xjiji for all i and j}. 

Thus A is the diagonal of P^ x P^, while D is a divisor of type (1, 1). 
In particular, if we let tti and 7T2 be the projections P^ x P^ — >• P^ 
and let if be a hyperplane in P^, then D is linearly equivalent to 
nlH + n;H. 

We also observe that A C D. It follows from [11] that 



Aa(P,v) < Xd{P,v) + 0^{1) for all P e ((P^ x P^) \ \D\){K). (6) 



Here 0^,(1) denotes an Mi<-bounded function (i.e., it is bounded by 
an Mft-constant) in the sense of Lang [8]. Note that our assumption 
on A and B ensures that {A, B) is not in the support of D. We also 
note that since D is an effective divisor, the local height \d ^^ bounded 
below by an M^-constant for all P not in the support of D. Hence 
evaluating (6) at P = {A, B) and summing over all nonarchimedean v 



and assume that A^ B. Then 




yields 



aa((A5),^) < (AD((A5),^)+a(i)) 




< Y Ai.((A5),^)+0(l) 



v^Mk 



^hD{{A,B))^0{l) 
^h^*H+.iH{{AB))+0{l) 
= /i(A) + /i(S) + 0(l). 



(7) 
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It remains to compute the local height relative to the diagonal. 
(In [11], Aa is called an arithmetic distance function.) Prom [11] and 
knowledge of the generators of the ideal defining A, we see that a rep- 
resentative local height function for A is given by 

\/\i{A,B),v) — min viaiBj — ajBi) — min viai) — min viBi). (8) 

' > 0<i<j<N ^ ^'^^^ 0<i<N ^ ^' 0<i<N ^ ' 

Next we observe that for any nonzero ideal C = (71, ... , 7^), we have 
^^logNC= 

(The 1/[A' : Q] in front comes from the way that we have normalized 
the valuations in order to simplify the formula for the absolute height.) 
Hence when we sum (8) over all nonarchimedean places of X , we obtain 

^ Aa((^, B),v) = T^^aog NS - log N21 - log N«). 

Substituting this into (7) yields the desired result. □ 

We conclude this section with an elementary result saying that every 
point in P^(i^r) has integral homogeneous coordinates that are almost 
relatively prime. 

Lemma 7. Let Kj Q be a number field and let Rk be the ring of integers 
of K. There is an integral ideal C — €{K) so that every P e F^{K) 
can be written using homogeneous coordinates 

P = [ao,ai, . . . ,q;jv] 

satisfying 

aQ, . . . , G Rk (^n-d {chq, . . . , ajv) | C 

Proof. Fix integral ideals ai,...,a/i that are representatives for the 
ideal class group of Rk- Given a point P e P^(i^), choose any homo- 
geneous coordinates P — [/Sq, . . . , Pn]- Multiplying the coordinates by 
a constant, we may assume that f3o, ■ ■ ■ , (3n ^ Rk- The ideal generated 
by /3o, . . . , /?Ar differs by a principal ideal from one of the representative 
ideals, say 

(7)(/3o, ---,Pn)^ <Xj for some 7 e X*. 
Then each 7/3j e aj C Rk, so if we set ctj = 7/3i, then 

P — [ao, . . . , QiAr] with ao, . . . , £ Rk and («0) - - - , ctiv) = 

Hence if we let € be the integral ideal € = 0102 • • • a/i, then C depends 
only on K, and for any point P we have shown how to find homogeneous 
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coordinates in Rk such tliat tlie ideal generated by the coordinates 
divides □ 

3. An ideal characterization of orbit size 

In this section we estimate the size of a certain ideal having the 
property that its prime divisors are the primes with irip < m. 

Proposition 8. With notation as in Section 1 and Theorem 3, there 
are a constant C = C {K, V, tp) and a finite set of exceptional primes 
S — S{K, V, (fi) so that for every m> 1 there exists a nonzero integral 
ideal Dim) satisfying the following two conditions: 

(i) For p ^ S we have rup < m if and only if p | S)(m). 

(ii) loglogNS)(m) < Cm. 

Remark 9. If y is projective and (p is finite of degree d > 2, then the 
following more precise version of (ii) holds: 

log log NX) (m) < (log d)m + Clog m. 

Proof. By definition, mp{(p, P) is the smallest value of m such that 
there exist r > 1 and s > satisfying 

r + s = m and (P) = (p^{P) (mod p). 

Notice that s is the length of the tail and r is the length of the cycle 
in the orbit O^^ (P mod p). 

We let € be the ideal described in Lemma 7. Then for each n > 
we can write 

^''\P) = [Ao{n),A,{n),...,AN{n)] 
with Ai(n) e Rk and such that the ideal 

2l(n) := (Ao(n), . . . ,A]\[{n)) divides the ideal €. 
It follows that for all primes p f € we have 
(^'■+*(P) = (p'{P) (mod p) 

«^=4> Ai{r + s)Aj{s) = Ai{s)Aj{r + s) (mod p) 

for aU < i < j < iV. 

Hence if we define ideals ®(r, s) by 

Q3(r, s) = {A{r + s)Aj{s) - A{s)Aj{r + s))^^.^.^^ 
and define D{m) to be the product 

S)(m) = Yl ®('^'«)' 

r>l,s>0 
r+s=m 
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then for all primes p \ € we have 

mp{ip,P)<m <^=^ p I D(m). 

Thus that 1D(m) has property (i). Further, the assumptfon that P has 
infinite 99-orbit tells us that 

(p''+'{P) ^(f'{P) for allr > 1 and s > 0, 

so D{m) ^ 0. 

It remains to estimate the norm of S)(m). We apply Proposition 6, 
which with our notation says that 

Using the fact that N2l(r + s) and N2l(s) are smaller than N£, which 
only depends on we find that 

1 



■ log NQ3(r, s) < h{ip'+'{P)) + h{ip\P)) + 0(1). 



[K: 

Next we apply Proposition 4 to estimate the heights, which gives 

logNQ5(r,s) < Cd'■+^ 

where C = C{K, V, (p, P) and d = d{V, ip) > 2. The key point is that 
neither C nor d depends on r or s. 

The ideal D{m) is a product of various Q5(r, s) ideals, so we obtain 

logND(m) = J2 logN23(r,s) < ^ Cd^+* < Cmrf™+^ < Cd^"". 

r>l,s>0 r>l,s>0 
r+s=m r+s=m 

Taking one more logarithm yields 

loglogNS)(m) -C m, 

where the implied constant is independent of m. □ 

An immediate corollary of Proposition 8 is a weak lower bound 
for nip. 

Corollary 10. With notation as in Section 1, there is a constant C — 
C{K, V, (f, P) so that 

mp(</7, P) > Clog log Np for all primes p. 

Proof. For each m > 1, let D{m) be the nonzero the ideal described in 
Proposition 8. Then for all but finitely many primes p we have 

p I D(mp) and log log ND(mp) < Cm^. 



10 



JOSEPH H. SILVERMAN 



It follows that ND{mp) > Np, which gives the desired result for p. 
Adjusting the constant to deal with the finitely many excluded primes 
completes the proof. □ 



4. An analytic estimate 

We now prove the key analytic estimate required for the proof of 
Theorem 3. This analytic result is a dynamical analog of a theorem of 
Romanoff [10], see also [2, 3, 7, 9]. The proof, especially insofar as we 
obtain an explicit dependence on s, follows the proof in [9]. 

Theorem 11. With notation as in Theorem 3, let X> 1. Then there 
is a constant C — C{K, V, ip, A) so that 

E ^^<4 foralls>0. (9) 

Proof. To ease notation, define functions g{t) and G{t) by 

g{t) = and G{t) = e-''\ (10) 

We use Abel summation to rewrite the series S{(p,P,X,s) in (9) as 
follows: 



Np • e*"*^ 



pGSpec Rx 

= J2 ^(Np)G'K) 



peSpccF/f 



= E( E ^(Np)G(m)) 

m>l peSpeciix 
mp=m 

= j]G(m)( J2 9m- E am] 

m>l pgSpccRif pgSpcc_RK 

mp<m mp<m—l 

= ^(G(m)-G(m + l)) J] ^(Np). (11) 

m>l 

The mean value theorem gives 



m>l peSpec Rk 



A-l„-s6»^ 



G{m) - G{m + 1) < sup -G'{0) = sup sXe^'^e 

m<d<m+l m<d<m+l 

— sXm'' 
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Substituting into (11) yields 



m>l pGSpecRif 
mp<m 



Si^,a,X,s)<j:sXm^-h--' J2 i^. (12) 



To deal with the inner sum, we use two results. The first. Proposi- 
tion 8, was proven earlier. The second is as follows. 

Lemma 12. LetK/Q be a number field. There are constants ci andc2, 
depending only on K , so that for all integral ideals D we have 
v-^logNp , , 

J]^-^<Cil0gl0gNS) + C2. 

Proof. This is a standard result. See for example [9, Corollary 2.3] for 
a derivation and an explicit value for ci. □ 

Using the two lemmas, we obtain the bound 

S{ip,a,\s) 

<Y^sXm'-^e-">' J: from (12), 

m>l peSpec ilif " 

mp<m 

— sXm^~^e~'^"^^ from Proposition 8 (i), 

m>l peSpec Rk 

p|S(m) 

< ^ sAm'*'~'^e~^"'^ (ci log log ND(m) + C2) from Lemma 12 

m>l 

< CsX ^ m^e"^'^^ from Proposition 8(ii). 

m>l 

It remains to deal with this last series. If A = 1, then we can explicitly 
evaluate the scries, but this is not possible for general vahics of A. (For 
example, if A = 2, then it is more-or-less a theta function). Instead we 
use the following elementary estimate. 

Lemma 13. Fix A > and /jl > 0. There is a constant C — C{X,n) 
so that 

00 

J2 m^'e-''^^ < Cs-^^+^^Z^ for alls>0. 

m=l 

Proof We estimate 

V m^e-"™' < / t^e-^'*' dt 
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r"00 



M'^e""^ dt letting u = s^^H. 

The integral converges and is independent of s. □ 
Applying Lemma 13 with n — X and substituting in above yields 

S{ip, a, A, s) < Ci{K, V, ip)sX ■ C2(A)s-(^+i)/^ = Cs{K, V, ip, X)s-'/\ 
This is the desired result. □ 

5. Proof of Theorem 3 

We now use the analytic estimate provided by Theorem 11 to prove 
our main density results. 

Proof of Theorem 3. (a) For any < 7 < 1 we let 

= {p e Spec(i?i^) : mp < (logNp)^}. 
Then for all s > we have 

— > — log Np from Theorem 11 with A = I/7, 

peSpeciJx IMP • e P 

^ ^ logNp 



Hence 



> T. — ^ , ^ M by the definition of V~^, 

_ ^ logNp 



S{V'y) = limsup(s — 1) ^ ^ by definition of upper density, 



— (Np)^ 



^-^ logNp u ^1 

= limsup s > ... , replacmg s by s + 1, 

< limsup Cs^"''' from (13), 

s-*0+ 

= since 7 < 1. 

Since the density is always nonnegative, this proves that S{V^) = 0. 
This is equivalent to Theorem 3 (a), which asserts that the complement 
of Vj has density 1. 

(b) The proof is similar. Let e > and define 

Ve^{pe Spec{RK) : mp < elog Np}. 
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Applying Theorem 11 with A = 1 and using the definition of Ve, we 
estimate 

y. logNp ^ logNp ^ log Np 



g - . Qsm^ - /l^ |\|h . Qsm^ - ^ |\|h . e^elogNp 



^ ^ logNp 
(Np)i+^^- 

Replacing s by s/e yields 

Ce logNp 
(Np)i+*' 

Hence 

5{Ve) = Umsups < Ce. 

It follows that the complement of Ve has lower density at least 1 — Ce. 

□ 



6. Conjectures and Experiments 

The density estimate provided by Theorem 3 is probably far from 
the truth. If the (/9-orbit of a point P in V^(Fp) is truly a "random map" 
from V^(Fp) to itself, then the expected orbit length nip should be on 
the order of ^J^V{¥p) ^ Npa'i™^. (See [4, 5] for statistical properties 
of orbits of random maps and [1] for the analysis of orbits of certain 
polynomial maps.) Thus we might expect the quantity 



1 log Np 



logX ^ rn^'^'^^ 

° Np<X '"'P 



(14) 



to be bounded as X — > oo. (Note that Snp<x ^^Np^ ~^ ^0 

We tested this guess numerically for the case of quadratic polyomial 
maps on Pq by computing the value of (14) for various polynomi- 
als (p{z) — z"^ + c and using values of X encompassing every 200th 
rational prime up to 20000. The results are listed in Table 1. 

The data in Table 1 is ambiguous; all five columns are increasing, 
but it is not clear whether they are approaching a finite bound. We 
extended the calculations for + 1 up to X = 50000, and the value 
of the sum (14) continues to grow as shown in Table 2. This suggests 
that rup grows somewhat slower than \/Np. 

Emboldened by these few experiments, we are led to make the fol- 
lowing plausible conjecture. 
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X 


-2 


z'-l 




z' + l 


z^ + 2 


1223 


1.4733 


1.6042 


1.4156 


1.3539 


1.3533 


2741 


1.7770 


1.6576 


1.5116 


1.4089 


1.4232 


4409 


1.9864 


1.6937 


1.5711 


1.5165 


1.4911 


6133 


2.1050 


1.6964 


1.6860 


1.6068 


1.5751 


7919 


2.1657 


1.7330 


1.8091 


1.6314 


1.5988 


9733 


2.2507 


1.7372 


1.8555 


1.6596 


1.6148 


11657 


2.2868 


1.7622 


1.8825 


1.7212 


1.6722 


13499 


2.3366 


1.7782 


1.9351 


1.7223 


1.7049 


15401 


2.3928 


1.8822 


1.9973 


1.7307 


1.7226 


17389 


2.4279 


1.9119 


2.0528 


1.7376 


1.7475 


19423 


2.4551 


1.9211 


2.0726 


1.7421 


1.7582 



Table 1. ^ Ep<x ^ various ^{z) 



X 


6133 


13499 


21383 


29443 


37813 


46447 


Avg 


1.6068 


1.7223 


1.7627 


1.7790 


1.8092 


1.8398 



TABLE 2. ^ Ep<X S ^(^) = + 1 



Conjecture 14. Let K/Q be a number field, let (f : be 
a morphism of degree d > 2, and let P G P^(A') be a point with 
infinite (p-orbit. Then for every e > the set 

{p:m,{if,P)>Np--'} 

is a set of density 1 . 

It is reasonable to pose a stronger question. Is there some k > so 
that the conjecture is true for the set of p satisfying 

N 

rrip > Np2 (logNp)-''? 

And in the other direction, we ask whether for some constant C > 0, 
the set of p satisfying rup > CNp^^^ is a set of density 0. 
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